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ABSTRACT 


The  visioplasticity  method  of  analysis,  employed  extensively 
in  investigations  related  to  metal  processing  problems,  is  applied  to 
an  analysis  of  a  rigid  wheel  traversing  a  soft  soil  where  flow  is  con¬ 
fined  to  two  dimensions.  Emphasis  is  placed  on  the  way  data  taken  from 
X-ray  photographs  are  reduced  in  order  to  obtain  the  velocity,  strain- 
rate  and  strain  patterns  for  this  flow  condition.  Consideration  is  giver 
to  modeling  the  dominant  constitutive  or  material  behavior  and  how  this 
is  incorporated  along  with  the  governing  continuum  field  equations  to 
provide  calculations  of  the  energy  dissipation  rate  and  the  stresses 
acting  on  the  wheel.  Experimental  data  obtained  by  McGill  University 
is  used  to  illustrate  the  feasibility  of  this  approach. 

This  method  of  analysis  appears  to  be  ~  very  promising  tech¬ 
nique  for  investigating  soft-soil  mobility  problems.  The  need  for  better 
experimental  controls,  particularly  in  soil  preparation  and  simultaneous 
flow-field  visualization,  and  possible  future  applications  of  the  technique 
are  indicated. 


FOREWORD 


The  research  reported  herein  was  performed  by  Cornell  Aeronautical 
Laboratory,  Inc.  (CAL),  and  is  a  segment  of  the  Soil  Mechanics  Tas>'  of  the 
Off-Road  Mobility  Research  (ORMR)  program  sponsored  by  OSD/ARPA  for 
which  CAL  is  the  prime  contractor  and  the  U.  S.  Army  Research  Office  - 
Durham  is  the  contracting  agency.  The  reported  study  is  part  of  a  broad 
research  program  designed  to  acquire  off-road  mobility  knowledge;  develop 
analysis,  prediction  and  decision  methodologies;  and  organize  knowledge 
and  methods  to  facilitate  use  by  military  planners,  vehicle  designers  and 
field  personnel  concerned  with  off-road  mobility.  Under  these  broad  objec¬ 
tives,  soil  mechanics  research  is  directed  at  developing  techniques  for 
analyzing  off-road  mobility  problems  related  to  soft-soil  conditions.  The 
present  study  demonstrates  the  application  of  a  technique  used  in  metal 
processing  analysis,  called  the  viscoplasticity  method,  to  a  two-dimensional 
rigid  wheel  in  soft  clay.  The  study  was  performed  during  the  period  November 
1967  -  September  1968,  under  the  guidance  of  Mr.  George  Bartleti,  ORMR 
Program  Manager.  Mr.  A.  N.  Tedesco  was  ARPA  Technical  Monitor. 

Methods  for  data  reduction  and  computer  coding  were  developed  by 
Messrs.  John  M.  Grace  and  Seymour  Samet,  and  Mrs.  Miriam  J.  Watson. 

All  experimental  data  used  in  this  study  was  obtained  under  ORMR  sub¬ 
contract  by  the  McGill  University  Soil  Mechanics  Group  under  the  direction 
of  Prof.  Raymond  N.  Yong.  During  the  course  of  this  study.  Professor  Yong 
has  also  provided  many  helpful  ideas  concerning  the  development  of  the 
analytical  techniques. 
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1.  INTRODUCTION 


The  role  of  soft  soil  in  the  generation  of  tractive  forces  acting  on 

r  1  21 

off-road  vehicles  has  been  under  investigation  for  a  number  of  years.  *■  ’  ^ 

The  problem  is  generally  recognized  as  one  having  mixed  boundary  conditions 
with  large  inelastic  sinkage  where  significant  energy  is  dissipated  within  the 
soil  medium.  However,  surprisingly  little  effort  has  been  directed  toward 
making  measurements  within  the  soil,  and  most  research  effort  is  restricted 
to  measuring  only  external  parameters  such  as  drawbar  pull,  torque,  normal 
load  and  slip.  Furthermore,  the  soil  properties  are  usually  described  as 
being  either  cohesive  or  frictional,  with  little  consideration  given  to  the 
soil  state. 

Most  attempts  to  analyze  these  problems  have  avoided  the  mathe¬ 
matical  difficulties  by  making  some  far-reaching  assumptions  concerning 
the  soil  behavior  under  off-road  loading  conditions.  Although  these  assump¬ 
tions  may  make  the  problem  mathematically  tractable,  they  inevitably  lead 
to  engineering  models  which  do  not  satisfy  basic  continuum  mechanics  field 
equations.  In  this  sense,  these  engineering  models  do  not  explain  the  physical 
mechanism  responsible  for  the  generation  of  tractive  forces  on  the  vehicle. 

If  these  problems  are  approached  in  a  purely  analytical  manner,, 
it  soon  becomes  apparent  that  suitable  analytical  techniques  are  not  likely 
to  he  developed  in  a  reasonable  time  period.  As  a  result,  Cornell  Aeronautical 
Laboratory,  under  the  Soil  Mechanics  Task  of  the  Off-Road  Mobility  (ORMR) 
Program,  proposed  to  investigate  these  problems  on  a  somewhat  fundamental 
level  with  a  limited  effort,  guided  by  two  prime  considerations;  (a)  the 
solutions  must  be  based  on  continuum  mechanics  principles  with  recognition 
given  to  ihe  particulate  nature  of  the  material,  and  (b)  the  t  chniques  which 
are  developed  must  be  generally  applicable  to  all  types  of  off-road  loading 
conditions .  ^  ^  ^ 

Within  this  framework,  it  was  decided  that  the  probability  of  suc¬ 
cess  would  be  much  higher  if  combined  analytical,  empirical  techniques 
were  employed.  The  results  obtained  during  the  past  fifteen  years  with 
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the  visioplasticity  method  in  metal  processing  analysis  were  particularly 
impress!'  Furthermore,  it  appeared  that  this  technique,  along 

with  suitable  modification,  could  be  employed  in  investigations  of  typical 
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off-rc»d  loading  situations.  ^  ^ 

The  visioplasticity  method  assumes  that  the  velocity  field  can 
be  determined  experimentally,  and  that  the  stresses  can  be  obtained 
analytCically  from  these  measurements.  That  is,  the  strains  and  strain 
rates  are  calculated  directly  from  the  velocity  field  while  the  stress 
distribution  results  from  an  analytical  solution  uf  the  equilibrium  (or 
motion)  equations  and  the  constitutive  equation  relating  the  stress  and 
strain  for  the  material. 

This  technique  has  been  extensively  applied  in  metal  processing 
analysis  where  the  velocity  fields  are  determined  by  photographing  grids 
attached  to  the  surface  of  the  workpiece.  Photographs  are  usually  taken 
after  successive  incremental  displacements.  The  material  is  assumed 
to  behave  as  a  Levy*Mises,  rigid,  perfectly  plastic  solid. 

The  application  of  this  technique  to  soil  mechanics  problems 
requires  further  experimental  effort  in  three  areas:  (a)  measurement 
of  the  displacements,  (b)  preparation  of  the  soil  medium,  and  (c)  the 
determination  of  the  soil  mechanical  properties.  McGill  University  has 
successfully  developed  and  employed  the  pulsed  X-ray  technique  to  measure 
soil  displacements.  ^  ^  In  addition,  some  success  using  a  technique 

similar  to  the  visioplasticity  method  has  been  reported  at  Cambridge  University 
where  foundation-type  problems  were  investigated.  ^  Soil  preparation 

and  the  determination  of  the  mechanical  prop«:rties  requires  a  continuing 
effort  and  has  been  under  study  at  Cornell  University.  ‘ 

The  present  soil  mechanics  work  has  been  restricted  to  the  behavior 

of  nearly  saturated  clays  at  high  water  contents.  Because  such  materials 

generally  exhibit  little  internal  friction,  a  plasticity  theory  would  be  expected 

to  be  most  promising  with  this  type  of  soil.  Fortunately,  such  soils  are  of 
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prime  interest  in  soft-soil  mobility  problems,  '  ' 


The  purpose  of  the  present  study  was  to  determine  if  the  visio- 
plasticity  method  is  a  reasonable  way  to  analyze  the  typical  off-road 
loading  situation.  A  rigid  wheel  was  selected  for  the  study  because 
steady-state  behavior  could  be  assumed  and  McGill  University  was 
prepared  to  make  the  measurements  on  such  a  wheel.  It  was  initially 
recognized  that  full  utilization  of  this  technique  would  eventually  require 
the  instantaneous  visualization  of  the  entire  flow  field  and  a  better  speci¬ 
fication  of  the  soil  mechanical  properties.  However,  it  was  felt  that, 
rather  than  go  to  the  expense  of  constructing  elaborate  facilities,  it  might 
be  wise  to  use  existing  facilities  and  techniques  to  test  the  validity  of  this 
method.  In  this  sense,  the  results  of  this  initial  investigation  must  be 
viewed  in  a  qualitative  manner. 

If  it  is  possible  to  use  this  technique  to  analyze  these  problems, 
then  the  solutions  should  be  exact  in  the  engineering  sense.  Such  solutions 
can  serve  as  the  basis  for  verifica.tion  and  formulation  of  engineering  models. 
The  present  report  concerns  itself  vdth  an  illustration  of  how  the  visioplasticity 
method  can  be  applied  to  a  specialized  situation.  At  a  later  date,  the  results 
obtained  from  this  study  will  be  compared  with  those  of  other  specific  engi¬ 
neering  models. 
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2.  EXPERIMENTAL  MEASUREMENTS  AND  DATA  REDUCTION 
2.  I  Experimental  Measurements 

All  experimental  data  used  in  this  study  were  developed  by  McGill 
University  under  ORMR  subcontract  and  a  detailed  discussion  of  their 
experimental  facilities  is  available  in  Reference  11.  To  set  the  stage 
for  the  subsequent  analysis  and  discussion,  only  a  brief  description  of 
the  soil  preparation  process  and  the  type  of  experiment  are  presented. 

To  develop  the  soil  test  bed  with  lead  markers  in  the  wheel  mid¬ 
plane,  a  dry  kaolinite  is  placed  in  alternating  layers  with  a  prescribed 
amount  of  water  (approximately  50  percent  by  weight)  and  allov  .d  to  rest 
until  the  water  migration  is  equilibrated.  The  soil  test  bed  is  then  built 
up  from  this  material  in  layers  through  a  process  of  remolding  small 
quantities  and  tamping  them  into  position.  It  should  be  noted  that  the  soil 
is  built  up  so  that  the  plane  of  the  markers  is  horizontal.  The  lead  markers 
are  embedded  when  the  layering  is  one-half  the  required  depth.  When  the 
test  bed  preparation  is  completed,  the  bed  is  rotated  so  that  the  marker 
plane  is  vertical  and  in  the  midplane  of  the  wheel.  Walls  are  maintained 
along  the  sides  of  the  soil  so  as  to  prevent  any  side  flow. 

A  carriage  traveling  at  a  constant  velocity  carries  a  rigid,  13.  5- 
inch-diameter  wheel  having  a  width  of  3.  75  inches,  equal  to  that  of  the  soil 
test  bed.  In  addition  to  the  linear  velocity,  an  angular  velocity  is  imposed 
on  the  wheel.  Transducers  located  on  both  the  carriage  yolk  and  the  wheel 
axle  shaft  measure  the  drawbar  pull  and  torque  acting  on  the  wheel.  Sinkage 
is  recorded  by  measuring  the  relative  dista^ice  between  the  wheel  center  and 
side  rails. 

Before,  during,  and  after  a  test  run.  X-ray  photographs  are 
taken  for  a  6-  by  7-inch  area  (in  the  marker  plane).  Three  X-rays  are 
shot  during  the  test  with  each  shot  separated  by  the  time  required  for 
the  wheel  to  traverse  6  inches.  Thus,  in  a  coordi’-.ate  system  moving 
with  the  wheel,  these  three  X-ray  shots  map  out  the  soil  displacements 
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over  an  area  approximately  18  by  7  inches.  In  this  system,  the  X-ray 
measurements  are  equivalent  to  the  displacement  time  history  for  the 
soil  motions. 

It  is  readily  seen  that  the  marker  measurements  are  taken  in 
a  system  fixed  to  the  X-ray  unit  or  fixed  in  physical  space.  In  this  system, 
the  soil  motion  must  be  viewed  as  unsteady.  If,  however,  the  coordinates 
of  the  markers  are  viewed  relative  to  a  coordinate  system  attached  to  the 
wheel,  then  it  seems  reasonable  to  assume  a  steady -state  motion  for  con¬ 
stant  wheel  velocities.  Since  steady-state  motion  has  been  assumed  in  the 
analysis,  the  measurements  are  transformed  to  coordinates  attached  to  the 
wheel  axle.  This  simple  transformation,  along  with  data-s moothing  tech¬ 
niques,  is  shown  in  Appendix  I. 

2.  2  Calculation  of  the  Velocity  Field 

The  marker  positions  relative  to  the  wheel  are  shown  in  Figure  1. 

In  this  analysis,  we  assume  that  the  markers  shown  in  any  row  are  deformed 
from  an  initial,  horizontal,  straight  line  of  markers.  Consistent  with  the 
steady-state  assumption,  any  row  of  markers  can  be  taken  as  either  a 
streamline  or  a  path  line  for  a  material  particle.  Movement  along  a  row 
from  right  to  left  is  equivalent  to  following  the  motion  or  time  history  of 
a  particle.  Since  the  time  increme.its  between  markers  may  be  determined 
from  the  distance  between  these  markers  in  the  initial  or  undeformed  con¬ 
figuration  and  the  wheel  axle  velocity,  the  relative  time  history  for  any 
marker  is  known. 

The  results  shown  in  Figure  1  rc  ,'resent  some  smoothing  and 
the  precise  procedure  is  given  in  Appendix  T.  Nevertheless,  these  coor¬ 
dinate  points  are  used  to  determine  the  velocity  field.  As  shown,  a 
triangular  array  of  markers  was  employed  in  the  marker  plane.  Shown 
in  Figure  2  is  a  cutout  indicating  both  marker  points  and  velocity  points. 

If  we  consider  the  point  (i,j  )  and  its  neighboring  points,  then  the  x  and  y 
components  of  the  velocity  vector,  u  and  i'' ,  respe''tively,  are  defined 


5 


I 


at  the  point  according  to 


At: 


At- 


(2.  1) 


The  time  increment,  At^  ,  is  determined  from  the  initial  coordinate 
distances ;  i.  e.  , 


At  -  r 

‘  K 


(2.2) 


where  the  superscript  <,  indicates  initial  positions  (from  first  X-ray  shot) 
and  /  is  the  wheel  axle  velocity. 

In  this  manner,  the  data  are  tabulated  in  a  checkerboard  pattern 
where  particle  positions,  defined  at  marker  points,  alternate  with  velocity 
component  points  defined  at  midpoints  between  successive  markers  in  the 
same  row.  The  velocities  are  now  used  to  determine  strain  rates  and 
strains  at  the  marker  points. 

2.  3  Calculation  of  Strain-Rate  Fields 

Since  uniform  strain  conditions  are  assumed  throughout  the  soil 
thickness  and  no  side  flow  is  permitted,  the  only  nonzero  components  of 
the  strain-rate  tensor  are  6^  ,  6^  ,  and  6^^  .  In  terms  of  the  velocity 
components,  u  and  i/-,  these  components  are 
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(2.  3) 


Thus,  if  the  first  derivatives  of  the  velocity  components  are  known  at  the 
marker  points,  then  the  strain-rate  tensor  is  also  known. 

The  triangular  array  was  selected  primarily  because  of  these 
derivatives.  For  example,  with  a  central  difference  approximation, 
the  partial  derivatives  of  u  with  respect  to .«  and  y  at  the  point  (i^j ) 


% 


(see  Figure  2)  are  taken  as 


(2.4) 


Similar  approximations  could  be  made  for  the  derivatives  of  v'. 

Since  finite  strains  are  evident,  corrections  due  to  the  geometry 
changes  must  be  incorporated  into  the  above  approximations.  This  was 
done  by  fitting  surfaces  to  the  u  and  tr  velocity  components.  Some  smoothing 
is  also  employed  in  this  process  and  again  the  details  are  given  in  Appendix  I. 
Since  the  distance  between  marker  points  is  about  one  inch,  the  derivatives 
in  Equation  (2.4)  are  actually  determined  according  to 


y-;;) 

tf-  .  *0.S)  -  -  O.s) 


(2.5) 


where  the  u.  values  are  taken  from  the  surface -fit  function  about  the  point 
(c,j)  for  the  M  velocity  component.  It  is  clear  that,  with  this  scheme,  the 
velocity  values  arc  interpolated  near  the  velocity  points. 

For  the  analysis,  an  equivalent  strain  rate  (see  Section  3)  is  required. 
In  the  plane  strain  problem,  the  eqnivalent  strain  rate,  Z  ,  is  given  by 


*y 


(2.6) 


It  should  be  noted  that,  because  polynomials  up  to  second  order  are  used 
in  the  surface  fit,  this  scheme  is  equivalent  to  differentiation  of  the  surface 
and  evaluation  at  the  point  . 
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This  quantity  is  essentially  the  second  invariant  of  the  strain-rate  tensor 
(within  a  constant  multiple)and  is,  consequently,  a  scalar  function  inde¬ 
pendent  of  coordinate  transformations. 


2.  4  Calculation  of  Equivalent  Strain  Field 


In  order  to  incorporate  the  stress-strain  properties  into  the 
analysis,  it  is  necessary  to  determine  an  equivalent  strain  parameter. 
By  definition,  the  equivalent  strain,  ,  is 


(2.7) 


where  the  integration, with  respect  to  time, is  taken  along  the  particle  path. 
This  integration  is  readily  performed  along  the  path  lines  or  streamlines 
as  determined  from  Figure  1. 
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3.  MECHANICAL  PROPERTIES  OF  THE  SOIL  MEDIUM 

i'  I 

J  3.  1  Soil  Behavior  and  Loading  Conditions 

}  ^  In  a  system  of  continuum  mechanics  field  equations,  the  consti- 

I  I  tutive  equation  is  a  tensor  equation  which  describes  the  mechanical  prop¬ 

erties  of  the  medium.  Various  physical  behavior,  such  as  elastic,  plastic 
‘  I  or  viscous  properties,  can  be  modeled  by  appropriate  terms  in  the  con- 

■  stitutive  equation.  However,  since  precise  accounting  of  all  these  phe- 

>  I  nomena  leads  to  an  extremely  complicated  system  of  field  equations,  a 

realistic  application  of  continuum  mechanics  to  most  problems  requires 

■  that  only  the  more  dominant  aspects  of  the  material  properties  be  modeled. 

The  dominant  mechanical  properties  exhibited  by  materials  result 

■  from  the  material  state  and  the  loading  conditions.  Specifically,  fine  grain 

1’  "  soils  are  a  problem  to  military  vehicle  mobility  when  they  are  at  high  water 

•  f  14l 

H  content  and  near  full  saturation.  ^  •*  Furthermore,  the  loading  on  the  soil 

I  occurs  over  a  short  time  period  (measured  in  seconds),  The  maximum  or 

limiting  stress  conditions  are  reached  causing  considerable  soil  flow  and, 

I  upon  unloading,  there  is  little  elastic  strain  recovery.  As  a  result,  it 

seems  appropriate  to  initially  restrict  primary  attention  to  a  two -phase 

I  water  and  fine-grain  particle  material. 

Is  In  considering  a  two -phase  medium  such  as  a  saturated  clay,  the 

►  I  important  factor  to  consider  is  the  relative  load  velocities  of  the  water  and 

solid  phases  under  the  loading  conditions.  The  small  amount  of  drainage 
I  I  which  occurs  under  rapid  loadings  indicates  that  this  relative  velocity  is 

small.  Consequently,  for  off-road  loading  rates,  it  would  seem  appro - 

■  priate  to  consider  the  soil  as  a  single,  continuous,  incompressible  medium 

■  and  employ  total  rather  than  effective  stresses. 

■  Since  the  loading  condition  subjects  the  soil  to  large  strains  and 

■  limiting  stress  conditions,  our  initial  approach  relies  on  a  purely  plastic 

_  constitutive  equation.  Because  of  the  high  water  content,  the  limiting 

I  stress  may  tend  to  depend  upon  the  rate  of  loading.  Such  materials 

exhibiting  viscous  effects  will  require  a  more  refined  theory  relying 
I  upon  viscoplastic  constitutive  theory. 
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In  the  following,  simplifiec  rigid-plastic  and  rigid-viscoplastic 
constitutive  equations  are  discussed.  These  include  Levy-Mises  plas¬ 
ticity  and  the  Bingham  viscoplastic  relations.  Although  each  of  these 
relations  has  been  applied  primarily  in  other  areas,  the  governing  con¬ 
ditions  for  the  fine  grain  soils  indicate  that  these  constitutive  theories 
should,  with  suitable  modification,  provide  reasonable  first-order  approxi¬ 
mations  for  modeling  the  material  behavior.  A  detailed  derivation  of 
these  equations  can  be  found  in  other  sources  (e.  g.  ,  see  Prager,  Refer¬ 
ence  16). 

3.  2  Levy-Mises  Plasticity  Equations 

The  Levy-Mises  plasticity  equations  assume  that  (a)  the  flow  is 
incompressible,  (b)  there  is  no  recoverable  strain  upon  unloading,  and 
(c)  the  directions  of  the  principal  axes  of  the  stress  deviator  tensor  and 
those  of  the  strain-rate  deviator  tensor  coincide.  The  tensor  equation 
becomes 

{  O  F<k* 

tf;  .  »  < 

(3.1) 

where  •  is  the  strain-rate  deviator  tensor,  S;;  the  stress  deviator 
tensor,  A  a  scalar  porportionality  function,  F  the  yield  function,  and  k 
the  yield  stress  in  pure  shear. 

These  equations  must  be  augmented  by  a  suitable  yield  condition. 
For  the  Mises  yield  condition,  the  function is  equal  to  the  second  invariant, 
Ji  ,  of  the  stress  deviator  tensor.  Thus, 

^  (3.2) 

is  the  generalized  representation  of  the  Mises  yield  condition.  It  should 
be  noted  that  the  yield  condition  could  be  replaced  by  some  other  constraint 
such  as  the  Mohr-Coulomb  or  Tresca  yield  conditions. 

* 

See  Appendix  II  for  brief  review  of  the  tensor  notation  and  definitions. 
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For  a  rigid,  perfectly  plastic  theory,  k  in  Equation  (3.  Z)  is 
constant  but  for  a  work-hardening  material,  k  depends  upon  the  strain 
history.  There  exists  considerable  evidence  that,  with  saturated,  fine 
grain  clays,  some  yielding  occurs  with  initial  loading.  In  addition,  with 
continued  loading,  the  yield  stress  increases  in  a  manner  similar  to  a 
work -hardening  material.  Consequently,  in  the  representation  of  this 
soil  behavior,  it  is  further  assumed  that  a  yield  surface  is  formed  upon 
initial  loading  and  that  this  surface  expands  with  continued  loading. 

As  a  result  of  this  latter  assumption  and  that  of  material  incom¬ 
pressibility  (i.  e.  ,  the  strain-rate  tensor  is  equal  to  the  strain-rate 
deviator  tensor).  Equation  (3.  1)  becomes 


with 

F  m  Jg  *  k,^ 


(3.3) 

(3.4) 


where  now  k  is  taken  as  a  function  of  the  strain  and  £•  is  the  strain-rate 
tensor.  The  scalar  function,  A  ,  in  Equation  (3.  3)  must  be  determined  by 
experimentation.  The  usual  procedure  is  to  use  the  tensor  scalar  product 
to  square  both  sides  of  Equation  (3.  3).  Thus, 


■ 


(3.5) 


where  repeated  subscripts  indicate  summation. 

Following  Hill  ^  we  may  define  an  equivalent  strain  rate,  t  , 
and  an  equivalent  stress,  <r,  as 


e 


i 


3 


a-  « 


S--  S  - 


(3.6) 


With  this  definition  of  ,  the  relationship  between  <r  and  is  r  ■ -^5^ 
and,  according  to  Equation  (3.4),  r  m  ilk. 
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As  a  result  of  these  definitions,  Equation  (3.  3)  becomes 


■‘j 


(3.7) 


It  is  necessary  to  determine  a  test  which  will  provide  an  equivalent 
stress  versus  equivalent  strain  curve.  For  metals,  this  relationship  is 
usually  determined  by  a  uniaxial  compression  test.  With  soft  soils,  however, 
this  kind  of  test  is  difficult  to  perform  and,  in  Appendix  III,  we  show  how 
the  results  of  a  typical  soil  mechanics  cylindrical  (triaxial)  test  can  be 
interpreted  as  an  equivalent  stress  versus  equivalent  strain  curve. 


A  typical  equivalent  stress  versus  equivalent  strain  curve  for  the 
remolded  kaolinite  used  in  the  soil  test  bed  is  shown  in  Figure  3,  Stress  > 
strain  results  were  obtained  by  McGill  University  for  cylindrical  samples 
under  undrained  conditions.  A  number  of  tests  were  conducted  and  some 
variation  was  experienced.  The  curve  shown  was  used  in  subsequent  analysis 
and  represents  a  medium  for  the  tests.  Since  the  analysis  requires  deformation 
rates,  the  "natural"  or  logarithmic  strain  has  been  used. 


For  the  analysis,  the  hyperbolic  strain  relation,  developed  for 
metals  by  Cox  and  used  extensively  by  Kondner  has  been  used 

to  fit  the  equivalent  stress  '^ersus  equivalent  strain  curve.  The  curve- 
fit  relation  for  the  results  shown  in  Figure  3  is 


<r 


_ £ _ 

O.  0^6  +  o.  278SS' 


(3.8) 


In  the  computer  calculations,  this  equation  was  used  to  relate  equivalent 
stress  to  equivalent  strain. 

3.  3  Bingham,  Viscoplastic  Equations 

The  constitutive  equation  for  the  incompressible  Bingham  solid 


is 


O  if  -ff,  *  k. 

.  5-  if  {Tt  >  k. 


(3.9) 
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where/cis  a  kinematic  viscosity  coefficient,  the  second  invariant  of 
the  stress  deviator  tensor,  k  the  static  shear  yield  stress,  the  strain- 
rate  tensor  and  the  stress  deviator  tensor.  Again,  if  yielding  is  assumed 
to  occur  on  initial  loading,  then  Equation  (3.  9)  becomes 


i  .  ..  - r 

a-  *  S/u.  Z 


(3.10) 


when  equivalent  stress  and  strain  rate  are  employed. 

It  should  be  noted  that,  for  no  viscous  effects  (i.  e.  ,  for  ^*0), 
Equation  (3.  10)  reduces  to  the  Levy-Mises  plasticity  relation.  Equation 
(3.6).  Consequently,  the  Bingham  relation,  as  given  here,  may  be  taken 
as  a  first-order  extension  into  the  dynamic  or  viscous  range  of  the  static 
plasticity  theory. 

There  are  a  number  of  other  possible  material  representations 

for  viscoplastic  solids.  A  clear  discussion  of  the  constitutive  theory 

for  these  materials  is  given  by  Perasyna  ^  The  selection  of  a 

suitable  material  representation  will  require  extensive  testing  of  the 

material  at  various  loading  rates.  Nevertheless,  in  view  of  favorable 

f  2i  1 

results  reported  concerning  mud-drilling  investigations  the 

immediate  extension  of  the  plastic  results  to  a  viscoplastic  theory 
should  undoubtedly  employ  the  generalization  of  the  linear  Binghanri  solid. 

In  this  report,  no  calculations  are  performed  using  a  viscoplastic 
theory.  It  appears,  however,  that  the  visioplasticity  method  could  be 
readily  extended  to  include  this  kind  of  dynamic  behavior. 
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4.  CALCULATION  OF  LOAD  VARIABLES 

Once  the  kinematic  variables  are  determined  experimentally  and 
the  soil  mechanical  properties  known,  the  visioplasticity  theory  may  be 
used  to  determine  the  loads  (forces  and  stresses)  acting  on  the  wheel.  In 
the  following,  consideration  is  given  to  both  the  determination  of  the  energy 
dissipation  rate  within  the  soil  medium  and  to  the  stresses  acting  in  the 
soil  >  wheel  interface. 


4,  1  Determination  of  Energy  Dissipation  Rate 

A*^, cording  to  the  limit  analysis  theorems  of  Drucker,  et,  al. 
the  following  inequality  holds  for  any  kinematically  admissible  velocity 
field 


£  <  J  a-6 


kavds 


(4.  1) 


where  £  is  the  rate  of  external  energy  input,  a-  the  equivalent  stress,  L  the 
equivalent  strain  rate  for  the  elemental  volume  dv,  k  the  shear  strength,  and 
AV  the  velocity  discontinuity  along  the  surface  9  .  With  the  visioplasticity 
method,  the  velocity  field  is  measured  and  is,  consequently,  the  actual 
velocity  field.  As  a  result,  the  equal  sign  may  be  taken  in  Equation  (4.  1). 
Furthermore,  with  the  rigid-wheel,  soft-soil  problems,  the  only  observable 
velocity  jump  occurs  along  the  soil  -  wheel  interface. 

For  the  rigid  wheel  case,  let 


(4.  2) 


where  the  rate  of  internal  energy  dissipation.  ^  ,  is  given  by 


r  £  dv' 


and  the  surface  dissipation. 


by 


(4.  3) 


(4.  4) 


This  would  require  sorrie  modification  if  a  viscoplastic  theory  w-ere  employed. 
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with  £  the  soil  shear  strength  and  du'  the  velocity  jump  between  the  wheel 
and  soil  surfaces.  Both  of  these  energy  dissipation  components  are  readily 
determined  once  the  stress -strain  and  kinematic  properties  are  known. 

The  energy  dissipation  calculations  for  three  <^ests  are  shewn  in 
Figure  4  and  Table  I.  Also  shown  are  the  energy  input,  as  determined  from 
the  measured  external  loads, and  velocities.  Since  the  linear  velocities 
do  not  vary  over  a  large  range  (see  Table  I),  the  energy  dissipation  is 
shown  as  a  function  of  slip.  However,  for  large  variations  in  the  linear 
velocities,  some  other  time-dependent  variable  would  probably  be  more 
appropriate. 


Significantly,  the  results  indicate  that  energy  is  dissipated  in 
two  ways  because  of  the  stress  and  strain-rate  fields  in  the  soil  medium 
and  as  a  result  of  a  velocity  discontinuity  between  the  soil  and  the  wheel 
surface.  In  addition,  the  energy  dissjpated  in  the  soil  medium  appears 
to  show  some  variation  with  slip  (at  almost  constant  linear  velocity)  while 
that  in  the  interface  is  very  dependent  upon  the  slip.  This  latter  phenomenon 
requires  further  investigation. 


For  the  wheel  having  a  constant  velocity,  the  energy  dissipation 
may  be  interpreted  in  terms  of  an  equivalent  drag  force  acting  upon  the  wheel. 
For  example,  if  we  divide  by  a  standard  length  per  unit  time  (i.  e.  ,  the  wheel 
linear  velocity),  then 


£ 


(4.  5) 


may  be  taken  as  an  equivalent  drag  force  acting  upon  the  wheel. 


It  is  desirable  to  have  the  actual  stresses  acting  in  the  soil  -  wheel 
interface.  In  the  following,  we  show  how  these  stresses  may  be  calculated 
for  the  rigid  wheel.  In  the  analysis,  we  permit  velocity  discontinuities 
across  viscous  layers  of  high  energy  dissipation  but  no  stress  disconti¬ 
nuities  are  permitted. 


4.  2  Calculation  of  Soil  -  Wheel  Interface  Stresses 

Within  the  soil  -  wheel  interface,  the  stresses  and 

are  determined  by  the  followiiig  equations  (see  Appendix  IV  for  the  derivation) 
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(4.6) 

where,  theoretically,  ,  tg  ,  Yro  ^  ^re  strain  rates  determined  from 
the  velocity  fields,  while  5^is  the  equivalent  stress  taken  from  the  stress - 
strain  curve.  Equation  (3.  8). 

From  the  kinematic  measurements,  ^  andF-are 

known  at  discrete  points  throughout  the  soil  medium.  For  the  polar  coor¬ 
dinate  system  used  in  the  derivation  of  Equations  (4.  6),  the  Cartesian 
strain-rate  components  are  readily  transformed  at  each  discrete  point 
to  the  polar  form  by 


cos  26 

stn  26 

'4' 

-silo  26 

cos2P 

- 

while  the  quantities  <£  and  a- are  invariants  and  remain  unchanged  with  the 
change  in  coordinates. 

The  same  general  surface  fitting  routine  which  was  used  in  the 
determination  of  the  kinematic  variables  is  again  employed  tc  interpolate 
the  desired  values  on  the  soil  -  wheel  interface.  For  this  fit,  triangular 
arrays  of  fifteen  points  are  used  to  successively  interpolate  the  desired 
values  at  about  one-inch  inteivals  in  the  soil  -  wheel  interface. 

The  stress  results  for  the  test  having  the  lowest  slip  in  Table  I 
are  shown  in  Figure  5.  The  routine  is  currently  being  run  for  other  tests, 
and  it  is  anticipated  that  they  will  be  reported  in  a  separate  report.  As 
can  be  observed  from  Equations  (4,6),  we  would  expect  the  results  for 
the  shear  stress,  ^re  •  to  probably  be  somewhat  better  since  integration 
is  not  involved.  The  determination  of  the  radial  stress,  ,  is  more 
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complicated  since  two  angles,  6o  *  must  be  determined  in  order 

to  perform  the  integration.  It  is  difficult  to  estimate  these  angles  using 
the  X-ray  measurements.  Furthermore,  the  results  for  cr}.  are 
dependent  upon  the  value  taken  for  Bg, .  These  kinematic  quantities  in  the 
soil  -  wheel  interface  might  be  better  determined  by  augmenting  the  X-ray 
measurements  with  photographs  of  surface  grids  through  glass  walls. 


Some  comments  concerning  the  calculated  stress  distribution 
are  in  order.  The  radial  stress  appears  to  reach  a  maximum  value  near 
the  point  where  the  soil  enters  the  interface.  This  is  in  contradiction  to 
the  plate  -  sinkage  models,  where  the  maximum  radial  stress  is  obtained 
at  the  point  of  maximum  sinkage.  The  shear  stress  takes  on  its  maximum 
value  near  the  bottom  of  the  wheel.  As  the  soil  enters  the  interface  region, 
this  maximum  is  reached  gradually  but,  once  reached,  it  tends  to  remain 
constant. 

[?-3] 

The  results  might  be  compared  to  those  of  Onafeko  and  Reece 
where  stress  measurements  were  made  on  a  rather  large  wheel  (49-inch 
diameter).  These  results  tend  to  exhibit  some  of  the  trends  measured;  e.  g.  , 
maximum  radial  stress  is  built  up  quickly  in  the  contact  region.  It  should 
be  observed,  however,  that  a  comparison  is  not  necessarily  appropriate 
since  the  Onafeko-Reece  results  are  for  a  much  larger  scale  and,  also, 
their  problem  has  three-dimensional  flow  contrasted  to  the  two-dimensional 
constants  considered  here. 


To  compare  the  stresses  with  the  measured  loads,  the  following 
equatior.s  were  used; 


7--  (6.7S) 


W  -  (6.7  S)  (3. 7S)  \^J  (<ri.  sir>  B  -  cos  9)  ct^ 

D  --{(o  75)(3. 7S)\  f  (  stn  9-0-^  cos  9)  d9  1 


(4.  8) 
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where  T  is  torque,  \f\f  is  the  normal  load,  and  D  is  the  drawbar  pull.  The 
calculated  values  were  106  inch-pounds,  50  pounds  and  zero  (actually 
0.  25  pound),  respectively,  for  the  torque,  load  and  drawbar  pull.  These 
correspond  to  measured  values  of  113  inch.pounds,  74  pounds,  and  a 
negative  3  pounds,  respectively. 

It  is  interesting  that  all  of  the  calculated  values  are  somewhat 
lower  than  those  measured.  In  addition,  the  energy  dissipation  for  this 
test  was  lower  than  that  determined  from  the  measurements.  This  possibly 
means  that  the  stress -strain  relation  should  be  different  than  that  used  in 
the  analysis.  It  would  seem  that  this  supports  the  assertion  that  better 
determination  of  the  soil  properties  are  needed. 
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5. 


FUTURE  RESEARCH 


The  present  results  indicate  that  the  visioplasticity  method  can 
be  used  to  determine  either  the  energy  dissipation  rats  or  the  stresses 
within  the  soi'  medium.  However,  some  irr»provement  in  both  the  soil 
preparation  process  and  the  measurement  technique  is  necessary.  These 
problems  have  been  considered  and,  based  upon  the  results  of  the  Cornell 
University  study  various  ways  of  preparing  the  soil  medium  appear 

feasible.  The  problem  of  measurement  (namely,  the  instantaneous  visu¬ 
alization  of  the  flow  field)  may  be  obtained  with  commercially  available 
pulse  X-ray  equipment,  supplemented  by  surface  grid  photographs. 

The  visioplasticity  method  should  be  extended  to  other,  more 
important  loading  conditions.  Such  extensions  will  require  considerable 
effort  in  the  development  of  analytical  methods.  This  will  undoubtedly 
require  the  application  of  numerical  techniques  to  solve  the  governing 
equations. 

5.  1  Soil  Prepa-ation  and  Mechanical  Propfties 

At  present,  the  soil  behavior  in  mobility  research  is  assumed  to 
be  essentially  the  same  as  in  foundation  engineering  problems.  Because 
the  limiting  shear  stresses  exhibited  by  soils  used  for  foundations  are 
usually  much  higher  than  those  exhibited  by  soft  soil  common  to  vehicle 
impedance,  the  application  of  traditional  soil  mechanics  requires  con¬ 
siderable  extrapolation.  Recent  results  obtained  at  Cornell  University  ^  ^  ^ 
indicate  that  the  low  strength  behavior  of  highly  saturated  clays  is  con¬ 
siderably  different  than  anticipated  by  traditional  theory.  These  results 
indicate  that,  when  kaolinite  is  consolidated  to  a  low  stress,  it  exhibits 
reasonable  total  stress -strain  characteristics. 

When  considering  the  behavior  of  soft  soils,  attention  must 
necessarily  be  focused  on  fine  grain  soils  at  high  water  content.  These 
soils  tend  to  be  progressively  weakened  by  repetitive  loading.  Such  field 
behavior  evidently  results  from  considerable  soil  flow  or  shear  distortion. 
On  the  other  hand,  when  soils  are  normally  prepared  in  soil  bins,  the 
preparation  process  includes  combinations  of  soil  remolding,  tillage,  etc. 
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When  prepared  in  this  manner,  repeated  loading  will  tend  to  compress  the 
material  resulting  in  a  higher  shear  strength.  It  is  imperative  that  soils 
be  prepared  in  a  carefully  controlled  manner  so  that  they  exhibit  low  strength 
characteristics  and  virtually  incompressible  flow. 

It  is  recommended  that  some  effort  be  directed  toward  preparing 
test  beds  from  a  slurry  using  an  anisotropic  overconsolidation  process 
and,  furthermore,  that  the  mechanical  properties  be  expressed  in  terms 
of  total  stresses  rather  than  effective  stresses.  Although  considerable 
effort  would  be  required  to  prepare  this  kind  of  soil  bed,  the  mechanical 
properties  would  be  reasonably  well  defined  using  traditional  soil  mechanics 
equipment. 

Some  limited  success  might  be  obtained  through  the  use  of  artificial 
soils  where  the  fluid  phase  is  replaced  by  an  oil  or  some  other  medium. 

Such  materials  are  more  stable  with  time  and  might  be  more  easily  repro¬ 
duced  than  clay  -  water  systems.  This  type  of  material,  however,  should 
be  expected  to  provide  results  over  only  a  very  limited  range  of  soil  behavior. 
Even  though  artificial  soils  might  be  developed  so  as  to  simulate  the  required 
stress -strain  properties,  it  is  doubtful  that  they  could  simulate  the  strength 
changes  with  repeated  loadings. 

More  important,  clays  tend  to  take  on  water  Ly  swelling.  This 
builds  up  pressures  in  the  soil  which  causes  the  soil  to  be  left  in  a  state 
of  overconsolidation  upon  drying.  Such  loading  histories  are  perhaps 
most  responsible  for  the  state  of  fine  grain  clays  when  found  under  natural 
field  conditions.  It  is  highly  unlikely  that  artificial  soils  could  ever  simu¬ 
late  this  more  fundamental  behavior.  Thus,  it  is  reasonable  to  judiciously 
use  some  type  of  artificial  soil  in  laboratory  simulations  but  the  broader 
aspects  of  the  problem  will  require  experimentation  with  clay  -  water 
systems. 


Note,  effective  stresses  are  required  to  determine  the  degree  of  con¬ 
solidation  but  the  total  stress  concept  should  suffice  during  the  vehicle 
loading  simulation. 
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Some  additional  experimentation  is  required  to  determine  the 
stress -strain  properties  for  the  model  study  materials.  These  tests, 
employing  total  stresses,  should  be  performed  with  loading  rates  com¬ 
parable  to  those  observed  in  the  vehicle  loading  simulation.  In  the  rigid 
wheel  tests,  they  were  nominally  about  0.  25-inch  Anch/second.  Con¬ 
ventional  undrained  triaxial  (cylindrical)  tests  performed  near  this  rate 
in  both  the  confined  and  unconfined  conditions  should  be  satisfactory. 

If  the  soil  is  found  to  be  rate  dependent,  then  some  consideration  must 
be  given  to  a  viscoplastic  representation  of  the  material. 

5.  2  Types  of  Off-Road  Loading  Experiments 

Current  methods  of  visualizing  the  soil  flow  pattern  must  be 
improved.  The  present  method  employs  a  300-kv  X-ray  system  which 
permits  instantaneous  visualization  of  an  approximately  6-  by  7 -inch 
area.  Recent  studies  at  CAL  indicate  that  this  would  be  greatly  enhanced 
with  a  newer  600 -kv  system.  This  system  would  permit  the  visualization 
of  a  circular  area  approximated  by  30  inches  in  diameter. 

It  is  recommended  that  this  system  be  augmented  by  taking  photo¬ 
graphs  of  grids  placed  on  the  soil  adjacent  to  a  transparent,  lubricated 

wall.  The  McGill  University  group  has  observed  with  simultaneous  X-ray 
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and  photographic  measurements  that  the  wall  effects  may  be  important.  '■ 
Nevertheless,  the  problem,  particularly  in  the  soil  -  vehicle  interface, 
requires  some  high  speed  measurements  which  probably  can  only  be 
obtained  with  photographic  techniques. 

The  visioplasticity  method  should  be  extended  to  other  more 
important  problems.  We  cite  here  two  categories:  (a)  so-called  running - 
gear  elements,  and  (b)  soil  measuring  instruments.  Wheels,  tracks, 
grousers,  etc.  are  included  in  the  former  while  indentors  and  soil-surface 
measuring  devices  make  up  the  latter.  If  the  kinematics  of  each  problem 
can  be  adequately  determined  and  the  soil  properties  specified,  then  the 
visioplasticity  method  could  serve  as  a  powerful  tool  for  investigating 
the  behavior  within  the  soil  and  relating  the  mechanical  properties  to 
the  external  loads. 
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Although,  in  principle,  the  method  could  be  utilized  with  three- 
dimensional  problems,  present  investigations  should  be  continued  pri¬ 
marily  with  two-dimensional  experiments.  These  should  include  a 
flexible  wheel  and  a  wheel  -  belt  arrangement  simulating  a  track.  Such 
investigations  would  require  additional  analytical  effort  in  order  to  deter¬ 
mine  the  stresses.  For  example,  in  the  rigid-wheel  analysis,  the  stress 
solution  is  simplified  since  it  is  readily  established  that  the  radial  stress 
component  in  the  soil  -  wheel  interface  depends  only  upon  the  angular 
coordinate.  This  simplification  would  not  be  possible  in  an  investigation 
of  a  flexible  wheel  where  a  numerical  scheme  would  probably  be  required. 

Present  measuring  instruments  essentially  attempt  to  relate 
measured  forces,  torques,  etc.  to  limiting  shear  stresses.  Because  the 
boundary  value  problems  have  not  been  solved,  these  devices  correlate 
the  soil  mechanical  properties  with  measured  loads  in  a  purely  empirical 
manner.  A  more  precise  investigation  should  shed  some  light  on  the  kind 
of  correlation  and  the  number  of  tests  required  to  draw  the  necessary  con¬ 
clusions  regarding  the  soil  state. 

Finally,  it  should  be  observed  that,  even  though  the  prime  interest 
is  currently  with  saturated,  high-water  content  clays,  eventually,  some 
consideration  will  have  to  be  given  to  other  soils.  Such  investigations  will 
again  be  most  dependent  upon  the  way  the  soil  is  modeled  in  the  constitutive 
equation. 
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OQ 


Table  I 

COMPARISON  OF  MEASURED  AND  CALCULATED  ENERGY  DISSIPATION  RATES 


SLIP 

PERCENT 

LINEAR  VELOCITY 
in. /tec 

ENERGY  DISSIPATION 
in.-Ib/tec 

RATE 

MEASURED 

CALCULATED 

4 

£ 

1 

E 

3 

3.1 

65 

49 

5 

54 

19 

3.2 

80 

54 

15 

69 

H6 

4.1 

148 

75 

86 

16) 
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APPENDIX  I 


COORDINATE  TRANSFORMATIONS  AND 
DATA -SMOOTHING  TECHNIQUES 


'’-1.  Coordinate  Transformations 

There  are  actually  two  transformations  required  in  order  to 
place  the  data  in  the  assumed  steady-state  coordinate  system;  i.  e.  ,  from 
the  film  plane  to  the  marker  plane  and  from  the  fixed  marker  coordinate 
system  to  the  moving  system. 

A  schematic  representation  of  the  experimental  setup  is  shown 
in  Figure  I-l.  Since  all  measurements  are  made  on  X-ray  film,  they 
represent  a  magnification  of  the  actual  distances  in  the  physical  plane 
of  the  markers.  If  we  let  (Xf  .yf)  represent  points  in  the  film  plane  and 
cor  responding  points  in  the  marker  plane,  then  it  is  readily  seen 
that  the  required  transformation  is 


where 


(I-l) 


with  and  the  distances,  respectively,  from  the  X-ray  optical  center 
to  the  marker  plane  and  film  plane. 


This  transformation  is  followed  by  a  single  translation  from  the 
fixed  coordinate  system  of  the  marker  plane  to  that  moving  with  the  wheel 
denoted  by(';t,yl.  As  is  shown  in  Figure  I-l,  the  coordinate  distances  from 
the  wheel  axle  to  the  optical  center  are  given  by  Since  three 

X-ray  shots  were  taken  during  the  test,  Of{t^)  takes  on  three  values  whose 
successive  difference  is  6  inches,  while  remains  constant.  At  any 

given  time,  the  following  transformation  applies: 


H  --  (*,) 


iO 


(1-2) 


1-2. 


Data  Smoothing  of  Flow  Pattern 

It  is  assumed  that  there  are  two  kinds  of  inaccuracies  in  the  data. 
The  analysis  requires  that  all  points  initially  lie  on  a  straight,  horizontal 
line  and  the  first  X-ray  shot  indicates  a  slight  departure.  Furthermore, 
inaccuracies  result  from  the  basic  problems  of  measurement  and  control, 
particularly  in  the  matching  of  the  results  from  one  shot  to  the  next. 

From  the  measurements  of  the  first  X-ray  film  (undisturbed 
markers),  the  y  coordii  ate  values  are  averaged.  Then,  for  each  marker, 
the  deviation  from  this  mean  is  determined.  This  deviation  is  then  added 
to  the  y  coordinate  value  of  the  same  marker  in  successive  X-ray  film. 

It  is  believed  that  this  adjustment  represents  a  first-order  correction  which 
tends  to  account  for  small  departures  from  a  purely  horizontal  alignment 
in  the  initial  configuration. 

Since  the  relative  times,  ,  are  known  for  all  markers  in  the 
X-ray  photographs  taken  during  the  test,  the  coordinate  values  of  each 
marker  are  expressed  as  Thus,  the  measured  values  at 

each  marker  are  expressed  as  a  function  of  the  time  f.  For  a  marker 
point  in  any  row,  a  second-order  polynomial  is  fitted  to  the  x  values  and 
also  the  y  values  using  -t  as  the  independent  variable.  Five  points,  centered 
about  the  marker  point,  are  employed  with  additional  constraints  in  the  least 
squares  sense.  The  corrected  value  is  obtained  by  evaluating  the  curve  fit 
function  at  the  known  value  of  .  This  scheme  is  applied  to  all  data  points 
on  the  three  X-ray  photographs  taken  during  the  test. 

1-3.  Velocity  Interpolation  Scheme 

Since  the  velocity  components  are  not  available  on  a  regular  grid, 
some  interpolation  is  required.  Furthermore,  it  was  found  that  some 
data  smoothing  was  required  because  of  the  matching  problem  and  also 
to  permit  more  flexibility  in  the  interpolation  subroutine.  A  general 
surfac e -fitting  subroutine  was  developed. 
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For  example,  if  we  wish  to  fit  a  function  fPfz.y))  about  the  point 
fz'.yV.  let 


P(x.i^)  •  (*•*")  *  Am( •4-y*)  *  As^x-r*)  ^ 

a(x-x*)  .  (1-3) 

Thus,  at  least  six  data  points  are  required  and,  if  more  are  given,  the 
normal  equations  may  be  developed  using  traditional  least  squares  con¬ 
straints.  The  subroutine  is  organized  so  that  any  number  of  points  may 
be  specified.  In  the  surface  fitting  for  the  velocity  components,  where 
derivatives  are  desired  at  the  marker  point  (c,J),  the  twelve  nearest 
neighboring  velocity  points  were  specified.  This  selection  of  points 
is  modified  slightly  on  the  boundary  and  adjacent  rows  where  one-sided 
derivatives  must  be  calculated. 
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APPENDIX  II 


TENSOR  NOTATION  AND  DEFINITIONS 


A  more  complete  description  of  the  properties  of  the  stress,  strain 
and  strain-rate  tensors  can  be  found  in  most  continuum  mechanics  textbooks. 
In  the  following,  however,  we  present  a  very  brief  review  of  some  of  the 
more  important  properties  and  definitions  for  the  stress  and  strain  tensors. 

Stress  and  strain  are  second-order  tensors  and,  for  example,  the 
stress  tensor,  oj^-,has  nine  components  since  i  and  J  each  take  on  values 
one  through  three.  These  indices  correspond  to  spacial  coordinates  such 
as  X,  y  and  ^  .  With  the  indicial  notation,  repeated  indices  indicate  sum¬ 
mation;  e.  g.  , 


<^kk  ~  °~3s  *  °~*  *■  ^  *  , 


(II-l) 


The  Kronecker  delta,  ,  is  defined  as  zero  when  L  and  j  differ 
and  unity  when  l  and  y  have  the  same  value.  Thus, 


■  ■  a  \ 


O  if  i 


t  if  i  J  . 


(II-2) 


The  following  definitions  are  made  for  the  stress  tensor,  o-^  j  , 


and  the  strain  tensor,  6—  : 


(a)  hydrostatic  or  mean  stress,  p, 

(b)  stress  deviator  tensor,  , 

iS  =  0~  *•  pS;- 


(II-3) 


(II-4) 


(c)  mean  strain,  e, 


(II-5) 
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(11-6) 


(d)  strain  deviator  tensor,  e 


tj  • 

e  ■  ■  -  ~  c6 ■  ■ 

tj  ‘j  •’j 


The  invariants  of  a  tensor  are  determined  from  its  characteristic 
equation.  For  example,  with  the  stress  tensor,  cr-  ,  the  characteristic 
equation  is  the  determinant  equated  to  zero, 


(T- 


(II-7) 


which  is  equivalent  to  the  cubic  equation 

*  I ,  er‘ I  ^  r  *■  I  ^ 


(II-8) 


The  coefficients  J^,  and  /j  are,  respectively,  the  first,  second  and  third 
invariants  of  the  stress  tensor.  The  principal  stresses,  trj  ,  oj ,  and  tr^  are 
the  roots  of  the  cubic  equation. 
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APPENDIX  III 


CYUNDRICAL  SOIL  MECHANICS  TEST 


In  the  cylindrical  or  so-called  triaxial  soil  mechanics  test, 
cylindrical  soil  specimens  are  subjected  to  controlled  longitudinal 
deformations  and  lateral  pressures.  For  the  undrained  test  on  satur¬ 
ated  clay  soils,  the  behavior  is  almost  incompressible.  In  order  to 
satisfy  the  requirements  of  the  visioplasticity,  it  is  necessary  to  inter¬ 
pret  the  results  in  terms  of  equivalent  stress  and  equivalent  strain 
parameters. 

A  cylindrical  coordinate  system  (r.9,^  )is  selected  such  that  the 
}  coordinate  corresponds  to  the  longitudinal  axis  of  the  cylinder.  Throughout 
the  test,  it  is  assumed  that  <r0  ,  O'^  and  are  the  only  nonzero 

stresses  and  strains,  respectively.  It  is  also  postulated  that 


% 


becomes 


(III-6) 


-  lia- 

>  *  ZW  > 


By  definition,  the  equivalent  stress  is 

or,  substituting  Equations  (III-l)  in  Eqioation  (III-7), 

<F  -  (a^  -  jt;.). 


(III-7) 


(III-8) 


As  a  consequence. 


(III-9) 


Thus,  during  the  test,  there  is  a  one-to-one  correspondence  between  the 
strain  increment  and  the  finite  strain,  and  the  curve  (er^  -^)  versus  6y  can 
be  taken  as  the  required  relation  between -^and 
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APPENDIX  IV 


CALCULATION  OF  STRESSES  IN  THE  SOIL  -  WHEEL  INTERFACE 

Consider  the  polar  coordinate  system  attached  to  the  rigid  wheel 
as  shown  in  Figure  IV- 1,  In  the  system  strain  rates,  6^  ,  60  >  and  , 
and  the  equivalent  strain  rate,  £  ,  are  determined  from  the  measured 
velocity  field.  The  equivalent  stress,  is  determined  from  the  stress - 
strain  curve.  A  solution,  in  the  visioplastic  sense,  means  that  the  stresses, 

<75  and  ,  must  be  determined  as  functions  of  thi‘se  measured  variables. 

The  plasticity  relations  for  an  incompressible,  Levy-Mis«*s  theory 

are 

‘  '2  ^  P) 

3  S  , 

eg  *•  j  ^  {cTo  *-p) 

(iv-D 

where 

”  "J  *  a-g  *  (T'y) 

The  hydrostatic  pressure,  p  ,  may  be  eliminated,  yielding 

‘=n- 

=  5  f-  yr0.  (IV-i) 

The  engineering  shear  strain  rate,  Yrt  ^  used  here;  i.  e.  ,  '^^r0  • 
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gives 


The  equilibrium  equations  in  polar  coordinates  are 


3<r^  O~r-<^0 

dr  r  d&  r 


*  L  ^  Hl 

r  r" 


(IVO) 


Substitution  of  Equations  (IV-2)  into  the  first  equilibrium  equation 


rr.  2  t\r  ,  ■  ..  /^/c7=.  \ 

r  '  ~  J  *  2  ^a(t  ^'■7] . 


(IV-4) 


Integrate  to  obtain 


<7-,  •  -  Jj/-P[f  (^r  7^  (1V.5) 

where  ^  is  the  wheel  radius  and  the  integration  function,  ^(3),  is  the  value 
of  C7^  at  r*  Since  the  stresses  are  required  at  the  interface,  the  problem 
is  reduced  to  determining  K(6), 

To  determine  K(6),  differentiate  the  first  of  Equations  (IV-2)  with 
respect  to^,  use  Equation  (IV-2)  to  eliminate  (^  ,  then  substitute  this  equation 
along  with  the  second  plasticity  equation  into  the  second  equilibrium.  Equa¬ 
tions  (lV-3),  This  results  in 


2K(0) 


(IV-6) 


If  r-i?  ,  then  the  integral  in  the  above  equation  is  zero.  Solving  for  K(6),  yields 
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The  function,  tCfiJ  ,  is  obtained  through  integration  with  respect  tod  and 
observing  that  <7^  s  ^7  at  r  .  /?  and  0  •  ;  i.  e.  , 


fCC0)  = 


f  cr  . 

7 1 


de 


i 


(IV-8) 


As  a  result  of  Equations  (IV-5)  and  (IV-2),  the  stress  in  the  soil 
wheel  interface  are  given  by 


rrjf 


K(e  ) 


/•--e 


r  X  f 


R  . 


Figure  H-l.  COORDINATE  SYSTEM  USED  FOR  STRESS  ANALYSIS 
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